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I. INTRODUCTION 



CD ■ 

Field theoretical models in external electromagnetic fields are of great interest and have attracted a lot of attention 
! They are relevant to the study of many physical systems whose properties depend on the effects produced by 

■ external fields, especially when such effects are non-perturbative in nature. A particular example of such a situation 



r—i 

in 



JlfJ (sec also |llj-|l3|] for some earlier studies) which 
|16| and in studies of the early Universe jl7) . 



is the so-called magnetic catalysis of chiral symmetry breaking 
CS| . has many potential applications in condensed matter physics [[14 

It is well known that a sufficiently strong attractive interaction between massless fermions can result in a chiral 
symmetry breaking condensate and a dynamically generated fermion mass. In the Nambu-Jona-Lasinio (NJL) model, 
for example, this occurs when the interaction strength is at a particular large critical value. It is also now known 
that the presence of an external magnetic field can be a strong catalyst for this effect, leading in some cases to the 
generation of a dynamical mass for fermions even by very weak attractive interactions. For example, in an external 
Q\ I magnetic field, the critical value of the NJL coupling is reduced to zero. 

It was suggested in Q that this magnetic catalysis of chiral symmetry breaking is a rather universal phenomenon 
and that its main features are model independent. The arguments were based on the observation that the effect is 
primarily due to the dynamics of the fermions from the lowest Landau level which are gapless and should therefore 
q; dominate the behavior of the system at long wavelengths. 

In this paper we will develop this idea further. We shall investigate the low-energy dynamics of chirally symmetric 
' models with a short range interaction in an external magnetic field. The idea is to construct an effective action 
| for the low energy degrees of freedom. In strong external magnetic fields, these are the modes of charged particles 
in the lowest Landau level. In the absence of interactions, their effective dynamics is described by a field theory 
whose dimensionality is two less than the spacetime dimension. We shall take the hypothesis that, at least for strong 
magnetic fields and short ranged interactions, the dynamics of chiral symmetry breaking can be understood in the 
lower dimensional effective theory. 

As an illustration, we shall argue that this assumption is very natural in 2+1 dimensions. There, the modes of the 
lowest Landau level are non-propagating, and in the absence of interactions, the ground state is highly degenerate 
and contains both chirally symmetric and chirally non-symmetric states. Interactions, even with infinitesimally small 
coupling constants, will split this degeneracy, and whether chiral symmetry breaking occurs or not depends on whether 
the interaction favors the chirally non-symmetric or chirally symmetric ground states. One by-product of our analysis 
will be the fact that the dynamical generation of parity violating masses seems to be suppressed by magnetic fields. 

In the more complicated case of (3+l)-dimensional gauge theory, we shall argue that the infrared dynamics is 
governed by an effective (l+l)-dimcnsional Gross-Neveu like theory with an infinite number of flavors of fermions. 
We derive the most general low energy effective action which is compatible with the symmetries of the four dimensional 
theory. We then study the resulting model using the renormalization group. We will present a plausible argument 
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for our main conclusion: that the presence of chiral symmetry breaking depends crucially on the sign (but not the 
magnitude) of one particular moment of the coupling constants which we denote by g v . If it is of the appropriate 
sign, its infrared renormalization group flow is to strong coupling and chiral symmetry is broken. With the other sign 
it flows to zero coupling and chiral symmetry remains unbroken. 



II. (2+l)-DIMENSIONAL PARADIGM 



Consider a system of relativistic (2+l)-dimensional fermions with the Dirac Hamiltonian 

„ _ ( ia-D 
H °-( -ia-D 



(1) 



where D = V + ieA(x) is the covariant derivative and a = (a , a 2 ) are the first two Pauli matrices. This Hamiltonian 
is obtained in the continuum limit of some tight-binding lattice models relevant to condensed matter physics — on a 
square lattice with half of a quantum of magnetic flux through each plaquette (see |18 and references therein) or on 
a honeycomb lattice without a magnetic field JlG 
the rf-wave state of high Tc superconductors ||14| 



.It has also been discussed in the context of the continuum limit of 



A. Symmetries 

Hq commutes with any linear combination of the matrices 

which generate the Lie algebra of a U(2) flavor symmetry which we will call chiral symmetry. Possible mass terms 
which can be added to H must be matrices which anti-commute with it. The basic matrix which anti-commutes 
with Hq and which commutes with all of the flavor generators is 

and a mass term which one could add to Hq which would preserve the flavor symmetry is j3m. 

Generally, in 2+1 dimensions, fermion mass terms violate parity symmetry. A parity transformation which is a 
symmetry of the massless Dirac equation 

H tP e (x) =E^ e {x) (4) 
and which commutes with the flavor symmetry generators T a is 

M*\x 2 )^K _1 2 ^je(-x\x 2 ). (5) 

The mass term j3m is odd under this transformation and it would therefore violate parity symmetry. More generally, 

all mass terms of the form (3 ^moI + m-Tj are odd under the parity defined in Eq. (j|). Moreover, there is no 

modification of the parity transformation which would make the mass 0m parity symmetric. 

It is possible to get a parity invariant mass by breaking the flavor symmetry. For this, we must define parity as a 
combination of the transformation (^) and a discrete flavor transformation, for example, 

M* 1 ,* 2 )- vj^feH 1 ,^)- (6) 

This is a symmetry of Hq and it also commutes with the mass term (3 {m 2 T 2 + m 3 T 3 ) . This mass term breaks the 
U(2) flavor symmetry to U(l) x U(l). 
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B. Zero modes 



This existence of the unitary matrix (3 which anti-commutes with Hq implies that Hq has a symmetric spectrum: 
if HipE = EipE, then Hf3ipE — —Ef3ipE and, with a particular choice of phase, ip-E — flipE- Thus, for each positive 
energy state, there is a negative energy state and the probability measure ip E (x)ipE(x) is a symmetric function of E. 

The interesting feature of Hq is the appearance of zero energy states in a background magnetic field. Consider the 
equation 

g ■ Duq = 0, (7) 
where the background electromagnetic field is static and has the property V x A = B(x). It has the solution 

uq(x) — exp [ie(f)(x) — ea 3 x{x)] v(x), (8) 

where 

-V 2 4>(x) =V ■ A{x), -V 2 x(i) = VxI(i)=B(i) (9) 

and 

a-Vv(x) = 0. (10) 

The second equation in Eqs. (j^) has the solution 

X(x)=-^ J d 2 y\n\x- V \B(y), (11) 



with asymptotic behavior 



lim X {x) = ~-\n\x\, x=^~ / d A yB{y). (12) 

|a;|^oo e Z7T ' 



(This is strictly correct if the magnetic field goes to zero sufficiently quickly at infinity for the total flux x to be finite. 
If it is not finite, there is an infinite number of zero modes and their density is given by x/V where V is the volume.) 
Asymptotically, 

fim u (x) = \x\ x °\{x). (13) 

In order to have a normalizable zero mode, the spinor v{x) should be an eigenvector of er 3 with eigenvalue — sgn(x) 
(— 1 if x > an d +1 if x < 0). It then must be a normalizable solution of the equations 

( d d \ 

a 3 v(x) = -sgn(x)v(x), f + isgn(x)^ j v(x) = 0. (14) 

Solutions of this equation are v(x) ~ [a; 1 — isgn(x)s 2 ] k for k = 0, 1, ... , [\x\] — 1 where [|x|] is the largest integer less 
than \x\ and the maximum value of k is the largest power allowed by normalizability of the wave function. Thus there 
are [|x|] normalizable zero modes. This fact depends only on the total magnetic flux and is independent of the profile 
of the magnetic field. The result, which we have found by explicit computation, is also a result of the Atiyah-Singer 
index theorem generalized to open spaces [^0| . Some other applications of this idea are discussed in |^l| . 
There are thus 2[|x|] zero mode wave functions of Hq. All of them obey the equation 

(3ipQ = -sgn(x)"0o, (15) 

[|x|] of them obey the equation 

PT 3 ^q = sgn( X )Vo, (16) 

and [|x|] of them obey 

PT^q = -sgn(x)Vo. (17) 

Of course, any linear combination of them is also a zero mode and could be chosen to be eigenvalues of any of the 
flavor symmetry breaking mass operators (3rh ■ T. All linear combinations would satisfy Eq. (051) . 
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C. Interpretation 



The existence of zero modes implies that the Dirac ground state of the second quantized system of fcrmions is 
degenerate. In the ground state, all negative energy states of the fermions are full and all positive energy states 
are empty. For overall charge neutrality, half of the fermion zero modes must be occupied — but the occupation of 
zero modes is otherwise unconstrained^] This means that the degeneracy of the ground state is Nq\/[(No/2)\] 2 where 
Nq = [\x\] is the number of zero mode wave-functions. 

In second quantization, the expectation value of any fermion bilinear operator, using the Dirac commutator for 
normal ordering, is [E| 

If d 2 x : *t(x)M*(x) :\ = -M ^ / d 2 xip E (x)Mil> E (x) - ^ J d 2 x^ E (x)M^ E (x) . (18) 

\ occupied unoccupied / 

Since, in the Dirac ground state, all positive energy states are unoccupied and negative energy states are occupied, 
and since (5 maps positive energy states onto negative energy states, if the matrix M commutes with /3, only the zero 
modes survive in this summation. 

First of all, this means that in a magnetic field the expectation value of the parity violating mass operator vanishes, 



J d 2 x : :\ = l^J d 2 x : & 0V :\ = 0, 



(19) 



regardless of the occupation of zero modes. This has the implication that external magnetic fields do not enhance the 
generation of parity violating mass terms. In fact, it indicates that if the external field is strong enough that the zero 
modes are well isolated, the expectation value of a parity breaking mass operator, (J : \f rT / 9 1 I r :) should always vanish. 
On the other hand, the parity invariant, chiral symmetry breaking mass (with M = /3T 3 ) has the expectation value 

J d 2 x : <b(x)TH(x) :\ = -isgn(x) I ]T J d 2 x^{x)T z M*) ~ E / d 2 x4{x)T^ (x) j . (20) 

\ occupied unoccupied / 

It is possible to choose ground states where this expectation value has any value between two bounds: 

- [|X|] < (/ d 2 x : *(x)T^(x) ^ < [\ X \]. (21) 

If we added an infinitesimal mass term m(3T 3 to the Hamiltonian Hq, the system would choose the ground state with 
minimal /3T 3 , that is the one where?] 

r d 2 x:*(x)T 3 *(x) =-[| X |]sgn(mx). (22) 

The observation that this is so for infinitesimal mass was made in However, this does not mean that chiral 
symmetry is necessarily broken. It only implies that the ground state is degenerate and there are chirally non- 
symmetric ground states which are degenerate with chirally symmetric ones. If there are additional weak interactions, 
of course, the correct ground state should be found by resolving the degeneracy using degenerate perturbation theory. 
In this case, all that is required to break chiral symmetry is an interaction which favors a chirally non-symmetric 
population of the zero modes. 



x Of course, occupying half of the zero modes is only possible when their number is even, which we shall assume. If it is odd, 
there are no neutral ground states. 

2 This is just the contribution of zero modes. If the total magnetic flux x is finite, the Dirac Hamiltonian has a continuum 
spectrum whose support begins at zero energy. Furthermore, there is a threshold density of continuum states whose contribution 
to bilinear densities such as Eq. (^) makes up the difference between [hd] and \x\ which is proportional to the spectral 
asymmetry of the Hamiltonian and is the correct result. For details, see [EOl. In the present paper, for simplicity we ignore 
the relatively small error which we make in neglecting the asymmetry of the continuum spectrum. When the magnetic flux is 
infinite, for example in the case of constant external field, [\x\] should be replaced by \eB\V/2TT where V is the spatial volume. 
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If the interaction favors the breaking of chiral symmetry, symmetry breaking occurs for even an infinitesimal value 
of its coupling constant. This is the reason why the critical value of the interaction can be at zero coupling in the 
presence of a magnetic field. In fact, in the absence of external magnetic fields a critical coupling necessary to break 
chiral symmetry is typically large. In the presence of the magnetic field it is reduced to zero. 

For example, it is well known that an interaction of the NJL model, 

L int -^(*T 3 *) 2 , (23) 

when added to the Hamiltonian will break the chiral symmetry if the coupling constant is attractive and is greater 
than a particular critical value. This is seen by analyzing solutions of the Schwinger-Dyson equation for a mass 
condensate. It was shown in B that, in the presence of a magnetic field, the critical coupling moves from some finite 
value to zero — even an infinitesimal coupling will break the chiral symmetry. That observation is consistent with 
our finding here. 

An interesting test of this idea would be to examine the effect on system with the interaction 

L int = |(^) 2 , (24) 

which can break parity by generating a parity violating fermion mass if k has the appropriate sign and is of large 
enough magnitude Ej|. Our analysis seems to suggest that a large magnetic field should in fact tend to increase the 
magnitude of the critical coupling. 



III. STRATEGY IN D=3+l 



In 3+1 dimensions, the situation is more complicated. The zero modes of the Dirac Hamiltonian in a background 
magnetic field still have some dynamics which is non-trivial and there are more possibilities for interactions. 

To understand the general strategy that we will take, begin with the Lagrangian which describes Dirac fermions 
interacting with an external electromagnetic field, 



So 



(25) 



where = dn + ieA^(x). For concreteness, we choose the chiral representation of the Dirac matrices, 



-i 

1 



a 1 
a 1 



(26) 



and we shall use the notation 



The Lagrangian has the form 



S a = J d A x i^l (d + #- D^j * L + (l> - a 



D) * 



(27) 



(28) 



This action has a global U(1)l x U(1)r chiral symmetry. In the following, we must keep in mind that when there 
are both background electric and magnetic fields, or if the gauge field is dynamical, this symmetry is reduced to the 
vector U(l) by the axial anomaly. For simplicity, in this paper we will assume that this is not the case. It would 
be straightforward to extend the present work to models with dynamical U(l) fields, by introducing more species of 
fermions so that there exist chiral symmetries which are unaffected by the axial anomaly. For now, we will consider 
the case of a constant non-dynamical background magnetic field, 



A Q =0, Ai = -—e lj3 x 



(29) 



In this background field, it is convenient to make use of the mixed spacetime-momentum representation as follows: 



5 



-^d 2 z [tft (w - fed 3 + X>) * L + (w + fca 3 - V) 9 R \ , (30) 
where, by definition, z = x 1 + ix 2 , d 2 z — dx 1 dx 2 and 

The spectrum of (fcer — X>) is well known. The equation 

(ka 3 -V)4>x{z,z) = \4> x (z,z) (32) 

has eigenvalues 



A = ±y/k 2 + 2n\eB\, (33) 

for all integers n = 0, 1, 2, . . . . When n > 1, these are dispersion relations of (l+l)-dimcnsional Dirac fermions with 
masses given by \j2n\eB\. [Of course, this mass gap is not Lorentz invariant from a four-dimensional point of view. 
However, it is Lorentz invariant in 1+1 dimensions. The reason for this is that there is a subgroup of the Lorentz 
group which survives in the background magnetic field — it is invariant under boosts in the direction of the field 
lines. Here this will mean that the effective theory for Landau levels must be invariant under the (l+l)-dimensional 
Lorentz transformations.] 

Zero modes of D are solutions of the equation 

V<f> m (z,z) = 0, (34) 
and are given by the infinite set of ortho-normal functions 

eB x(m+D/2 , eB v , 



= ^r(m + i) {-) exp {—**) [lj ' (35) 

with m = 0, 1, 2, . . . . Here, without loss of generality, we assumed that eB > 0. In the (l+l)-dimensional theory, the 
zero modes of T> correspond to massless fermions. The Lagrangian is (restoring the space dependence) 

/oo 
dx°dx 3 ( i ^ ) * d +^m ) + itf^-tf') + (massive modes), (36) 

m=0 

where d± = do =p ^3 and the new fields are defined as the coefficient functions in the expansions of and over 
the complete set of eigenstates: 

00 

*l = ^ ^m ) (a; ,2: 3 )'/ , m(z, z) + ^(massive modes), (37) 

771—0 

OO 

= ^ ^m\ x0 ,a; 3 )0 m (z,z) + (massive modes). (38) 



m— 



Note that the kinetic term in the action ( p6[ ) for the massless modes appears to have a U(N)r x U(N)l with 
N — ► 00 symmetry. This is just the unitary symmetry which mixes the different wave functions of the degenerate zero 
modes, and which is actually there for every Landau level. This effective symmetry is not preserved by interactions. 

We will analyze the possibility that interactions that are added to this field theory drive a spontaneous breaking 
of the U(1)ji x U(1)l chiral symmetry by generating a mass gap for the fermions in Eq. (|3^). We expect that this 
spontaneous symmetry breaking takes place at very long wavelengths. Our approach to the problem in the following 
sections can be summarized as follows: 

1. We consider the theory of four-dimensional fermions in a magnetic field as described above, with some inter- 
actions which should be local and respect (3+l)-dimensional Poincare and chiral symmetry but are otherwise 
unspecified. 

2. We then consider the effective field theory which is obtained by integrating out all of the massive modes of the 
fermions in Eq. <pw and momentum states of the massless modes above an ultraviolet cutoff. 
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3. We assume that the resulting effective Lagrangian is local. To guarantee locality, we generally have to assume 
that the interactions between fermions in four dimensions are short ranged]^ Furthermore, the ultraviolet cutoff 
for the effective theory should be less than the mass gap of the lightest massive mode, \eB\. 

4. We consider all relevant operators which could be added to the effective Lagrangian which are consistent with 
the symmetries of the theory. Since they should be relevant in the sense of two dimensional field theory, these 
are all possible four-fermion operators which are consistent with symmetry. There are an infinite number of 
such operators. 

5. We compute the beta function for the coupling constants of the relevant operators and look for infrared stable 
fixed points. These fixed points should govern the behavior of the very long-wavelength degrees of freedom of 
the theory. 

6. If the coupling constants flow to an infrared stable fixed point, then an infrared limit of the massless theory 
exists and there is no symmetry breaking. If, on the other hand, the coupling constant flow is not in the domain 
of attraction of any infrared fixed point, so that it flows to strong coupling in the infrared, we postulate that 
this implies the dynamical generation of a mass gap — and spontaneously broken chiral symmetry. We shall 
find examples of both kinds of behavior. 

An example of a local four dimensional interaction which preserves the U(1)r x U{1)l chiral symmetry is the NJL 
interaction p3], 

Si„t = ^J d4x [(**) 2 + (^75*) 2 ] ■ (39) 
A renormalizable version of this interaction would be one which is mediated by a massive scalar mesons. 

IV. GENERAL STRUCTURE OF THE LOW-ENERGY THEORY 

The constraints of (l+l)-dimensional Lorentz invariance and U(1)r x U(1)l chiral symmetry allow four-fermion 
coupling constants as in the effective theory, 

OC OO / \ 

£eff - E {t^WP + Wid-W) + E 90 ( £ 2 2 ) ( 4 °) 

n— ni,n2,mi,m2-0 

The coupling constants go ( ni 2 J obey further constraints from charge conjugation, parity and time reversal 



m± ni2 

(CPT) symmetry and the symmetry of the underlying theory under translations and rotations about the axis of the 
magnetic field which are summarized in Appendix [A]. A general solution of those constraints would yield the most 
general allowed structure of the low-energy effective action in Eq. ( pp| ) . 

Some particular solutions of those constraints are of interest. For example, the constraints allow the maximally 
symmetric solution with U(N)r x U(N)l (N — > oo) symmetry, 

9 ( Z 1 Z 2 ) = 9LRS ni , m2 8n 2 , mi , (41) 



mi m '2 

where gLR is real. Similarly, there is a solution with U(N)v (with N — > oo) symmetry. 



m n 2 
mi m,2 



— 9V^n\ ,712 ^Tn\ ,?7i2 5 (42) 



where gy is real again. 



li If there are long-ranged interactions which are mediated by massless fields, the correct procedure would be to retain the long 
wavelength modes of the massless fields in the effective Lagrangian. 
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In addition a large class of solutions could be found by the reduction of the interactions in the original (3- 
dimensional model to the lowest Landau level. After such a reduction of the NJL interaction (£59), we arrive at 



7ii n 2 
mi "*2 



G 



r(ni + mi + l)r(n 2 + m 2 + 1) 



2 (n 1 +mi+n 2 +TO 2 )/2 W p( ni + i)r( mi + l)T(n 2 + l)T(m 2 + 1)' 



(43) 



Since it descends from a Lorentz and chirally invar ia nt int eracti on in 3+1 dimensions, it must necessarily satisfy the 
constraints of symmetry [see Eqs. (A6), (|A7j ), ( |A10| ), ( All ) and ( |Al 2j )] . Unlike the first two solutions in Eqs. ( |4l| ) and 
(42|), this last one in Eq. (^) does not seem to have any extra symmetry in addition to the required U(1)r x U(1)l 
flavor symmetry. 

Now, we envisage having obtained the effective action ( f40| ) by integrating out all modes in the higher Landau levels, 
as well as all momentum modes of the fermions in the lowest Landau level which are above a certain cutoff. Of 
the many interactions that this procedure would produce, we have kept only the local four-fermion operators. This 
procedure is legitimate only if the ultraviolet cutoff of this model is lower than the lowest mass gap of the fields which 
have been eliminated, i.e. \eB\. By chiral symmetry and Lorentz invariance, the effective action cannot contain mass 
terms for the fermions. Furthermore, the only Lorentz invariant four-fermion operator is of the form given in Eq. pfj). 

The renormalization group procedure examines how the coupling constants in Eq. (40|) change as we further lower 
the cutoff to isolate the very long wavelength excitations. This information is encoded in the beta function. 

The j3 function for the general coupling constant in Eq. ( [4Cj| ) is computed in 2 + e dimensions to two loop order in 
Appendix |b[ The result is 



^ I Mi M 2 



2eg 



Ni N 2 
Mi M 2 



-E 



9 



Ni N 2 

fci Ni - N 2 + ki 



9 



Ni- N 2 + ki fci 
Mi M 2 



-y 



8tt 2 



fci ,fc 2 



+ g 



Mi Nx + Mi 

Ni k 2 

fci Ni + fci - fc 2 



fci k 2 

Mi Mi + fei-fc 2 



iVi + Mi-fei N 2 
fci M 2 



Ni + ki- k 2 N 2 

Mi M 2 + fci - fc 2 



M 2 



-fci 



fc 2 fci 
M 2 



Ni N 2 
fc 2 



- Mi 
M 2 



Mi + fci-fc 2 N 2 
N 2 -Mi + k 2 fci 



(44) 



(4tt) 2 



Nj. N 2 
Mi M 2 



E 

fci ,fc 2 



Ni fc 2 + Ni \ / fci - fc 2 fci 
fci fci - fc 2 J 9 { fc 2 + Ni Ni 



, fci k 2 - 
+ 9 [N 2 fci 



fci fc 2 + Mi 
Mi fci - fc 2 



9 



fci - fc 2 Mi 
fc 2 + Mi fci 



9 



M 2 k 2 - 
fci fci 



N 2 
-k 2 

M 2 
-fc 2 



fci - fc 2 N 2 
9 {k 2 + N 2 fci 



fci — fc 2 fci 
fc 2 + M 2 M 2 



Finally, before we proceed to the next section, we note that we may have to occasionally cut off the summation over 
modes in the first Landau level. We do this by summing modes up to some maximum number, n = 0, 1, . . . , N. Such 
a situation could be produced by considering the fermions in an external magnetic field with finite magnetic flux, or 
with uniform field and finite transverse area. In that case, N which can be thought of as the number of "flavors" 
which goes to infinity as \eB\Si 2 /2ir where Si 2 — > oo is the area of the two-dimensional perpendicular (with respect 
to the direction of the magnetic field) subspace. 



V. SYMMETRIC FIXED POINTS 



In Appendix |b[ we have computed the two-loop (3 function ( |i4| ) for the most general coupling constants, subject to 



the restriction in Eq. (A12). The next natural step would be locate infrared stable fixed points of this beta function. 
This problem is far more complicated than we can solve at present. Instead, in our analysis below, we shall restrict 
our attention to the beta function for coupling constants with close to maximal symmetry. 

For some scalar field theories, where there are many components of a scalar field and renormalizable interactions 
which couple them, it is known that the only infrared fixed points of the renormalization group flow of relevant 
couplings are those with maximal symmetry |23|,^4| . To our knowledge, there is no similar theorem for Gross-Neveu 
like models. On the other hand, we consider it plausible that similar arguments can be applied: in particular, if no 
maximally symmetric infrared stable fixed point exists, then there are no infrared stable fixed points at all. 



8 



Maximally symmetric couplings can easily be shown to be contained in other combinations of coupling constants. 
In Appendix |Cj we show that the NJL coupling contains the maximally symmetric ones. 

Consider the renormalization group flow in the special case when U (N)r x U(N)l and U{N)v couplings, as defined 
m Eqs. © and @, are the only non-zero couplings in the effective action. Then, from our general result in the 
previous section, we extract the following expressions for the (3 functions of interest: 



Plr = 2eg LR + ^g 2 v + (Ng 3 LR + 2g v g 2 LR + Ng 2 v g LR - g 3 v ) , (45) 



(3 V = 2eg v - ^-g 2 v + JL [ 2 Ng v - (N - 3)g 2 v g LR - (A - 2)g v g 2 LR ] . (46) 

Note that we have been forced to introduce a cutoff on N. In order to get a sensible result when taking the limit 
A — > oo, it is convenient to rescale the couplings as follows: gLR — > gir/VN and gv — > gv/N. As is clear from the 
definition, the (3 functions will also get rescaled accordingly, (3lr — > Pi r /V~N and /3y — > Pv/N. After taking this into 
account and performing the limit A — > oo, we arrive at 



Ar - 2e. 9ir + ^2 5? r , (47) 
(3 V = 2eg v - ^-g 2 v + ^fl flf r . (48) 



It must be emphasized that we have adopted the rescaling of the most symmetric gLR coupling by IfyN rather than 
1/N. Since the rescaling performed above led to the well-defined (3 functions, we conclude that the (l+l)-dimensional 
action in Eq. (Efl) with g^R of order ~ 1/yN describes a consistent and non-trivial interacting theory in the limit 

n -f oo. n 



When e = we can solve the renormalization group equations (47) and ( f48| ) explicitly, and the analytical solution 
reads 

9lr(t) = , 2 , (51) 

V 27T 2 

<?„ (t) = g " (0 ) . (52) 

This flow is presented graphically in Fig. [l] where the arrows show the flow direction toward ultraviolet. 

In the upper half-plane g v > 0, the simple analysis of the flow given by Eqs. ( pi] ) and ( [52] ) reveals infrared (with 
g v +00) and ultraviolet (with gi r — > ±00 and g v — > +00) Landau poles at 

= + (53) 



9v(0) 2 3 2( ) 

respectively. We argue that the strong infrared dynamics (with g v — * +00) in this half-plane of couplings is an 
indication of a mass generation and breaking of the chiral U(X) symmetry. Indeed, the generation of the fermion mass 
in the infrared region seems to be the only way one can avoid running into the problem of the physical Landau pole. 



4 As one can see, the rescaling of gLR by 1/N is also meaningful. The resulting f3 functions are 

Pir = legir, (49) 

1 2 

f3 v = 2eg v - —g v . (50) 

The corresponding theory is less interesting. Indeed, the expressions for the /3 functions in Eq. ( ^ ) and ( |5o| ) in the limit e = 
describe the situation when the gi r coupling does not run at all, while the g v experiences asymptotic freedom from the g v > 
side and infrared freedom from the g v < side. We shall see in a moment that this picture corresponds to a special case 
(gir = 0) of the flow described by Eqs. @ and @. 



9 



The generation of the fermion mass in the g v > half-plane, in its turn, is consistent with the expectation of the 
universality of the magnetic catalysis in a wide range of (3+l)-dimensional models (such as the NJL) with a short 
range interaction. Indeed, as we established, the low energy dynamics in such models is describ ed by t he (1 +1) 



dimen sional effective action in Eq. (|40|) with the coupling satisfying the set of restriction in Eqs. ( A10 ), ( All ) and 
( [A 12] ) . The generic coupling [say, like that in Eq. (|43| ) coming from the interaction of the lowest Landau level modes] 
which does not h ave a ny extra symmetry is still expected to have the U(N)r x U(N)l and U{N)y contribution [see 



Eqs. ( |C7aj ) and flC7b| )]. Then, if this U{N) V contribution is positive, g v > 0, it is going to drive the system to the 



generation of mass. 

Now let us study the flow in the lower half-plane g v < 0. As is easy to see, there is an infrared fixed point at 
(gir-,gv) = (0,0) and the ultraviolet Landau pole (with g v — > — oo while gi r is either fixed or approaches ±oo) at the 
following values of t: 



tuv = WWY if k(0)l < 4ir- (56) 

Since the infrared fixed point (0,0) corresponds to weakly coupled dynamics, there is apparently no mass generation 
in this half-plane of the coupling space. This is also in full agreement with our general expectation. Indeed, the 
negative values of g v correspond to the repulsion rather than attraction in the fermion- antifermion channel which is 
responsible for the generation of mass, breaking chiral U(l) symmetry. 

Our conclusion here is that whether this model breaks chiral symmetry or not is entirely dependent on the sign of 
the coupling constant g v . 

The reader might be puzzled by the fact that a U(l) chiral symmetry can be broken in an effectively two dimensional 
system. We emphasize here that this phenomenon is identical to that in the chiral Gross-Neveu model. Strictly 
speaking, chiral symmetry can only be broken in the large N limit. The finite N system should still be chirally 
symmetric, as the low dimensionality of the system would not allow for spontaneous breaking of a continuous symmetry. 
Indeed, if we consider the case of a large but finite value of N, an infrared stable fixed point appears for g v > at 
(gir,gv) — (—2ir,nN) which goes to infinity with N — > oo. This fixed point is at strong coupling, so a conclusion 
based on perturbation theory is speculative at best, but its appearance is consistent with the expectation that the 
massless limit of this model is well defined (albeit strongly coupled) when g v > and N is finite and chiral symmetry 
breaking need not take place. When N — > oo this fixed point moves to infinite coupling and chiral symmetry breaking 
is possible. 



VI. CONCLUSION 



In this paper we have shown that the effective low-energy dynamics of the U(l) chirally symmetric models with a 
short range interaction in a background magnetic field is described by a (l+l)-dimcnsional Gross-Neveu like model 
with an infinite number of flavors [see Eq. (|4Q)]. Here we established that different flavors come out as the represen- 
tation space of the magnetic translations in the original (3+l)-dimensional model. The number of flavors is infinite 
and proportional to the area of the two dimensional space perpendicular to the magnetic field. 

Ba sed exc lusiv ely on the arguments of symmetry, we established a set of conditions, given by Eqs. ( A6), (A7), ( [AlO] ), 
( All] ) and (A12), that the couplings of the effective theory have to satisfy. To show that they allow a non-trivial 
solution, we presented a few (out of infinitely many possible) examples of couplings that satisfy all the constrains. 
Among them, there are, in particular, the highly symmetric U(N)r x U(N)l and U(N)v (with N — > oo) couplings. 
These latter are of special interest because their renormalization group flow is self-contained and allows an analytical 
solution [see Eqs. (51) and (f>2|)]. 

At the level of the effective theory, we calculated the two-loop (3 function and analyzed the renormalization group 
flow in the two-dimensional subspace of the U(N)n x U(N)l and U(N)v couplings. The general result is argued to 
indicate the generation of the fermion mass in the g v > half-plane of couplings. In the other half-plane, the infrared 
dynamics is weakly coupled and there is no mass generation. This mass generation pattern is consistent with the 
earlier suggested conjecture of the universality of the so-called magnetic catalysis [|| in, at least, the models from the 
same universality class as the chiral U(l) NJL model. 
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APPENDIX A: CONSTRAINTS ON COUPLING CONSTANTS 

Let us clarify the origin of the action as well as the meaning of different flavors in Eq. (Eot) in terms of the original 



(3+l)-dimensional model defined in Eqs. fl2S| ) and (|39|) 



We start from the analysis of the space-time symmetries in the model. Notice that, due to the presence of the 
background field, the standard translations in the two-dimensional plane perpendicular to the magnetic field are not 
symmetries of the original model in Eq. (p8|). Nevertheless, there are other transformations, the so-called magnetic 
translations, which leave the action invariant. In contrast to the case of the ordinary translations, the two generators 
of the magnetic translations do not commute. In the symmetric gauge given in Eq. (|2j]), the explicit representation 
of the generators and their commutation relation read 

(Ala) 




(Alb) 

[X,Y] = i, (Ale) 

where we assume that eB > 0. It is easy to check that the operators X and Y commute with the Hamiltonian of our 
model. Instead of these X and Y, it is convenient to introduce the "creation" and "annihilation" operators 

a = , (A2a) 

* X-iY , A „, , 

at = (A2b) 

[a, at] = 1. (A2c) 

Having introduced these operators in the problem, we realize that the Fock space is spanned by the set of states \q, m) 
where the quantum number m — 0,1,2,... denotes the eigenvalue of a^a operator and the multi- index q represents 
all the other quantum numbers (say, the Landau level number, the fermion spin projection and the chirality). In the 
absence of any vacuum rearrangement (symmetry breaking), the above set of states (in coordinate representation) 
reads 



/I \ 1 1 

{x\m,n,a,x,P\\) = 



2ttZ (V2l) r 



^ ro " n) (f) exp (-f)^' (A3) 



where I = \j\feB is the magnetic length, z — x 1 + ix 2 , z — x 1 — ix 2 , X\\P\\ = x°pp — x 3 p3 and (f> atX is the spinor with 



a given spin projection a and chirality \. Note that the expression in Eq. (A3) is well defined even for the case when 



m < n (both numbers are positive) due to the Rodrigues formula for the generalized Laguerre polynomials, 

1 d n 

z m-n L {jn-n), y = e z — (e^Z™) . (A4) 

w nl dz n v ' y ' 

One has to remember also that all the modes in the lowest Landau level (n = 0) have the same projection of the spin, 
while the modes in the higher Landau levels (n > 1) have both projections. 

Since we are interested in the structure of the effective action, describing the infrared dynamics (po -C V eB), it 
is sufficient to take into account only those degrees of freedom that originate from the lowest Landau level modes 
(ri = 0). These modes freely propagate in the ( 1 + 1 )-dimensional parallel (x°,x 3 ) subspace, and are classified by 
the chirality and the eigenvalue of a^a [see Eq. flA3|)]. In what follows, we denote the effective degrees of freedom 
accordingly, ipj%' R \xu), where the superscript (L) or (R) denotes the states that result from the (3+l)-dimensional 
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states of definite chirality. In the parallel subspace (x°,x 3 ), ipm\%\\) and have the interpretation of the left 

and right moving along the x 3 axis modes, respectively. 

While restricting the kinetic term in the NJL model (|39|) to the lowest Landau level modes, 



_/ \ ^ (/> + ,+^\x\ l ) / z \ m ( zz\ 

^ii'^x: [wo cxp (-i^)' 

we check that 7^11^ — > 7|'p|i- After performing the integration over the perpendicular space coordinates in the original 
NJL action, we arrive at the effective model as in Eq. flip)) . Remarkably, while the NJL coupling G is dimensionful, 
the effective coupling in Eq. pp| ) is dimensionless, g ~ G/l 2 = G\eB\. 

From the derivation above, we see that the fields of different flavors in the effective (l-t-l)-dimensional theory (|4p| ) 
correspond to different eigenstates of a' a operator. This simple observation, as we show below, has far reaching 
consequences. 

Now, let us establish the allowed structure of coupling constants in Eq. ([l0|). The most general restriction on the 
couplings comes from the condition of reality of the action. This requires that 



ni n 2 \ _ I m 2 mi 
mi m 2 / \ nj ni 



(A6) 



Similarly, the invariance under the parity (x 3 — > —a; 3 ) leads to another restriction: 



I nx n 2 \ _ ( mi m 2 \ 
® \ mi m 2 J I nx n 2 J 

These two conditions are too general and so are not of great interest or of great power by themselves. It turns out, 
however, that there are other, more restrictive conditions. 

By recalling the origin of flavors in Eq. (^0|), we realize that the effective theory should enjoy some kind of flavor 
symmetry that results from the symmetry under magnetic translations of the original (3+l)-dimensional theory. This 
flavor symmetry, as will become clear in a moment, puts further restrictions on the allowed structure of the four-index 
coupling in Eq. (40). 

The infinitesimal transformations of the magnetic translations and the related rotation in the perpendicular plane 
are given by the following operators: 

Ux = l + iexX = 1 + ^p= (a + a f ) , (A8a) 
v2 

U 2 = 1 + ie 2 Y = 1 + i2= (a - a f ) , (A8b) 
v2 

U 12 = 1 + iei 2 a^a. (A8c) 

To determine the transformation properties of the fields of different flavors, we again recall that, by construction, 
these fields are the eigenstates of the a^a operator. Then, by doing a simple exercise, we find that the action of the 
creation and annihilation operators, cr and a, on the properly normalized fields should read 



a'ipn = Vn+ l^n+i) a>i>n = Vnip n -i- (A9) 

Making use of these properties, we check that the kinetic term in the effective action is invariant under the transfor- 
mations in Eqs. (A8a), ( A8b| ) and (A8c). In fact, if we had started with a more general, non-diagonal kinetic term 
in the effective action [which is not forbidden by the chiral £7(1) symmetry], the requirement of invariance under the 
magnetic translations would lead us back to the diagonal form as in Eq. (ppj) . 

Th e inv ariance of the four-fermion interaction g in Eq. (^0|) under the set of transformations in Eqs. (A8a), ( A8b| ) 
and (A8c) leads to 

— / ni - 1 n 2 \ , — p ( nx n 2 + l 

^{ mi m 2 J-^ Tl5 ( v m 1 m 2 

+V™Tg( ni , " 2 )-Vn^+Tg( Ul "«)=<>. (A10) 
v \ m i — 1 m 2 J V mi m 2 + 1 ' 

V^+Tg( ni + 1 n2 )~V^g( ni 712-1 

3 l mi m 2 J v I mi m 2 
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-\Jm x + Ig 



ni n 2 

mi m 2 



ni n 2 

mi + 1 771 2 



ni ri2 
mi m, 2 — 1 



= 0, 



= 0, unless ri\ — n 2 + mi — m 2 = 0. 



(All) 
(A12) 



APPENDIX B: CALCULATION OF THE TWO-LOOP /3-FUNCTION 

To derive the two- loop f3 function of the effective spinor theory in Eq. (f4(i|), we apply the method of Ref. [^4| that 
was used for the ip 4 scalar theory in d = 4 — s dimensions. In dimensional regularization (with D = 2 + 2s), our 
renormalized Lagrangian density reads 

00 00 , >. 

Cnjl = £ (zW^iftf^ + Z^^ ld .^) + ^ G ( Z Z ) 

71 — 711 ,fl2 ,™1 ,m2— 

OO / «. 

(V^-f^ + ^«L^) + E ^ ( £ m 2 ) ^?V£WVS2 + (• ■ 0. (Bl) 

n T Tin -m 1 1T> .-, O N ' 



= E . 

n—0 ni ,n2 ,mi ,m2=0 

where the coupling G includes the coupling constant renormalization, G — Z g g, and the ellipsis denotes the coun- 
tcrtcrms. We remind the reader that the fou r-index coupling g is non-zero only for n\ — n 2 + m\ — m 2 = 0. This 
means that the four-fermion term in Eq. (Bl) contains the sum only over three indices (say, n±, n 2 , and mi), while 
the fourth one {m 2 ) is superfluous. 

Before proceeding with the actual calculation of the (3 function, we need to specify how to handle the infrared 
divergences that appear in the calculation of the Feynman diagrams |24|^5|] . Such divergences usually come from the 
propagators of the massless fermions. If treated improperly, they could easily obscure the calculation of the relevant 
diagrams and eventually lead to a wrong result. To avoid the problem, in what follows, we modify the infrared region 
by changing the fcrmion propagators as follows: 

S^ L \p) = P —^^-^, (B2a) 

pZ pZ _ m z 

* (R) (p)=S-HT^r ( B2b ) 
p z p z — m 

This infrared regularization procedure respects all the symmetries of the model and does not change the ultraviolet 
region. 

Now, let us calculate the f3 function. First of all, we recall that the relation between the bare and the renormalized 
couplings, in the dimensional regularization (D = 2 + 2s), reads 

Ni N 2 \ fi- 2£ _ / Ni N 2 



90 { Mi M 2 J / (L) ( L) G { Mi M 2 J • ( B3 ) 

V A N t Z M X Z, N 2 Z M 2 

This is going to be used in the definition of the (3 function. In calculation, we impose the following renormalization 
conditions: 



9 A2R) 

d P J NN 



p=0 



d r (2i) 

d P+ NN 



r (4) 



( Ni 


N 2 \ 




\Mi 


M 2 ) 


= 9 

p=0 



= 1, (B4) 

Ni N 2 
Mi M 2 



(B5) 



Note that the introduction of the effective infrared cutoff m in the fermion propagators earlier allows us to use the 
renormalization point at p = 0. 

The Feynman diagrams of the relevant contributions to the two-point functions at two-loop order are given in 
Figs. ||a and ||b. By extracting the divergent (of order 1/e) terms of these two-loop corrections, we arrive at the 
equations 
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dp- 



-r 



(2R) 
NN 



dp+ NN 



p=0 



7 (R) 

J N 



7 {L) 
J N 



4r) 



(4tt) 2 £ 



E^ 

All ,k,2 



fci k 2 + N 
N ki- k 2 



~(0 



v ' y- r ( N k 2 + N 

(47r) 2 £ fe ^ \h ki-k 2 



G 



G 



kx-k-z N 

k 2 + N ki 

ki - k 2 ki 

k 2 + N N 



where, by definition, 
(Fig. 2a) -> 

(Fig. |b) 



4r) 



(4tt) 2 £ 



O(l) 



c9 



(4tt) 2 £ 



0(1) = - 



2D-4 
1 

(4-7r) 2 e 
1 d 



d D kd D q 



(k + p)-(q- k)+q- 



dp- J (2tt) 2D [(jfe + p) 2 - m 2 ] [[q - k) 2 - to 2 ] {q 2 - to 2 ) 



O(l). 



/./, 



2D-4 g p+ 
1 



d D kd D q 



(k+p) + (q - k)-q+ 



{2tt) 2D [(fc + p) 2 - to 2 ][( 9 - fc) 2 - m 2 ](g 2 - m 2 ) 



p=0 



(47T) 



0(1). 



Thus, we see that c' r ' = = — 1. 

In a similar way, the perturbative expansion for the four-point function reads 



r( 4 ) 

a + aie 



JVi iV 2 
Mi M 2 



G 



Mi M 2 



47T£ 



E 

fci 



G 



Ni N 2 

fc a iVi - iV 2 + fei 



G 



iVi - iV 2 + h h 
Mi M 2 



G 



iVi 



fci 



Mi iVi + Mi - ki 
b + b\e 



G 



Nx + Mt-h N 2 
fci M 2 



(47re) 2 



E 

fci ,fc 2 



G 



Ni N 2 
fci Ni - N 2 + fci 



G 



Ni- N 2 + fci 



fci 

Ni - N 2 + k 2 
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Mi M 2 



G 



Ni fci 

Mi Ni + Mi - fci 



G 
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iVi + Mi - fc 2 N 2 
fc 2 M 2 



'(4tt) 2 £ 



E 



ki,k 2 



G 



Ni fc 2 

fci iVi + fci - fc 2 



G 



Nt + h- fc 2 N 2 

Mi M 2 + fci - fc 2 



G 



+ G 
d + die 

' (47T£) 2 



fci fc 2 
Mi Mi + fci - fc 2 



G 



ATi N 2 - Mi + fc 2 
fc 2 M 2 



G 



Mi + fci - fc 2 N 2 
N 2 - Mi + fc 2 fci 



E 

k\ ,/C2 



G 



Ni fc 2 

fci JVi + fci - fc 2 



G 



Ni - N 2 + fci fci 
Mi M 2 



G 



M 2 + fci-fc 2 fci 
fc 2 M 2 



7V 2 

- iV 2 + fci 



Ni + ki- fc 2 

fc 2 JVi 



+ G 
+ G 
+ G 



fci fc 2 
Mi fci-fc 2 + Mi 

TVi Ni + fci - fc 2 
fci fc 2 

iVi + Mi - fc 2 iV 2 
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G 



G 



7V 2 - Ni + fci fci 
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Ni k 2 
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G 

G 



fci - fc 2 + Mi iV 2 - iVi + fci 
fc 2 M 2 

fc 2 fci 
Mi Mi + fc 2 - fci 

'fci fci + fc 2 - M 2 
fc 2 M 2 



(B6) 
(B7) 



(B8) 



(B9) 



(BIO) 



Note that the restriction ni — n 2 + mi — m 2 = is satisfied for each four-index coupling that appears here. The 
coefficients <Zj, 6j, Ci and di are defined by the following expressions: 



(Fig. 3) 



ao + aie 

47T£ 



0(£) 



,D-2 



d^fc 



(27T) B (fc 2 - TO 2 ) 2 
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and 



(Fig. 



(Fig. 



(Fig. 6) 



(Awe) 



(-0 



(47r) 2 e 



d + d\e 

(47T£) 2 



Here 7 sa 0.57 7 is t he Euler constant. 
From Eqs. (pl|) - (|B14D, we obtain 



47T£ 
,2D-i 



1 +e(l + 7) + eln 
d D kd D q 



Awy 2 
q 2 k 2 



0(e), 



(2w) 2D {q 2 -m 2 ) 2 {k 2 -m 2 ) 2 



(47T£) 
„'2 



l + 2e(l + 7) + 2eln 



0(e) = 



2D-4 



d D kd D q 



Awn 2 
k\q-{q-k) 



O(l), 



1 



2(Aw) 2 e 



(2w) 2D (k 2 - m 2 ) 2 (q 2 - m 2 )[(q - k) 2 - m 2 } 
+ 0(e), 



0(1) EE 



2D-4 



d D kd D q 



k 2 q-(k-q) + 



(Awe) 2 



(2w) 2D (k 2 ~ m 2 ) 2 (q 2 - m 2 )[(k - q) 2 - m 2 

0(1). 



ao = -1, 01 = -(1 + 7) -In 



6 = 1, h = 2(1 + 7) + 21n 



Co 



1 

2' 

do = ~\, 



47TM 2 ' 

9 

m 



After expressing the function G in terms of g, we arrive at 
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The (3 function is denned as follows 
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By making use of this definition, we calculate the two-loop (3 function, 
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where we already used the fact that a 2 , = bo = —2do which, by the way, is the necessary condition for the renor- 
malizability of the model. Note that the last term in Eq. (B18) app ears d ue to the renormalization of the two-point 
function. After taking into account the values of constants in Eqs. (B15a) - (B15d), we arrive at our final result for 
the P function: 
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APPENDIX C: NJL COUPLINGS CONTAIN MAXIMALLY SYMMETRIC ONES 



While having less symmetry, the coupling in Eq. ( (43| ) still could contain contributions of highly symmetric solutions 
as in Eqs. ( [51) and (|4^). In order to extract such contributions, we introduce the projection operators to the 
corresponding subspaces in the space of couplings, 



P( L «) [...]= lim 1^ 
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jy Tr (Lii) 



Tr( y ) 



p(V) 



lim 1^ 

N^oo 



N(N 2 — 1) 
iVTr (v) [■•■]- Tr (£fl) [■ ■ ■] 
N(N 2 - 1) ' 



where, by definition, 
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By applying the projection operators in Eqs. (CI) and 
contributions, 
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(CI) 
(C2) 

(C3) 
(C4) 

(C5) 
(C6) 



to the coupling in Eq. (f43|), we easily extract the symmetric 
2G 



N + 1 

2G 
N+l~ 



I {V) . 



(C7a) 
(C7b) 
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FIG. 1. The renormalization group flow in the {gir,g v ) plane. 
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FIG. 3. Diagrams contributing to r' 4 ^ at one loop which determine do and ai. 
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FIG. 6. Diagrams contributing to r' 4 ^ at two loops which determine do and d\. 



24 



